Evolution of spin-wave excitations in ferromagnetic metallic manganites 
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Neutron scattering results are presented for spin-wave excitations of three ferromagnetic metallic 
Ai- x A' x Mn03 manganites (where A and A' are rare- and alkaline-earth ions), which when combined 
with previous work elucidate the systematics of the interactions as a function of carrier concentration 
x, on-site disorder, and strength of the lattice distortion. The long wavelength spin dynamics show 
only a very weak dependence across the series. The ratio of fourth to first neighbor exchange 
(J4/J1) that controls the zone boundary magnon softening changes systematically with x, but does 
not depend on the other parameters. None of the prevailing models can account for these behaviors. 



Determining the evolution of the elementary magnetic 
excitations in Ai_ x A' x MnO^ (where A and A 1 are rare- 
and alkaline-earth ions respectively) is the first step in 
understanding the magnetic interactions in these doped 
perovskite manganites. According to the conventional 
double-exchange (DE) mechanism the motion of 
charge carriers in the metallic state of Ai- x A x MnC>3 es- 
tablishes a ferromagnetic (FM) interaction between spins 
on adjacent Mn 3+ and Mn 4+ sites. In the strong Hund- 
coupling limit, spin- wave excitations of a DE ferromagnet 
below the Curie temperature Tq can be described by a 
Heisenberg Hamiltonian with only the nearest neighbor 
exchange coupling 2]. At the long wavelength (small 
wavevector q), spin- wave stiffness D measures the aver- 
age kinetic energy of charge carriers and therefore should 
increase with increasing x 0, 0- While spin dynam- 
ics of some manganites initially studied appeared to fol- 
low these predictions 0, 0] , later measurements revealed 
anomalous zone boundary magnon softening deviating 
from the nearest neighbor Heisenber g H amiltonian for 
other materials with x ~ 0.3 |fl III ISl l2l llfll • Three classes 
of models have been proposed to explain the origin of 
such deviations. The first is based on the DE mechanism, 
considering the effect of the on-site Coulomb repulsion 3] 
or the conducting electron band (e g ) filling dependence 
of the DE and superexchange interactions [111 . The sec- 
ond suggests that magnon-phonon coupling |8lll2l| or the 
effects of disorder on the spin excitations of DE systems 
[l3j | is the origin for the zone boundary magnon soften- 
ing. Finally, quantum fluctuations of the planar (x 2 — y 2 )- 
type orbital associated with the A-type antiferromagnetic 
(AF) ordering may induce magnon softening as the pre- 
cursor of such AF order 01 ■ Although all these models 
appear to be reasonable in explaining the zone boundary 
magnon softening near x = 0.3, the lack of complete spin- 
wave dispersion data for Ai_ x A' x Mn03 with x < 0.3 and 



x > 0.4 means that one cannot test the doping depen- 
dence of different mechanisms and, therefore, the origin 
of the magnon softening is still unsettled. 

Very recently, Endoh et al. measured spin- 

wave excitations in the FM phase of Smo.ssSro^gMnOa 
(SSM045) and found that the dispersion can be de- 
scribed phenomenologically by the Heisenberg model 
with the nearest neighbor (Ji) and fourth-nearest (J4) 
neighbor exchange coupling (Fig. 1). By compar- 
ing the Ji/ Ji ratios, which measure the magnitude 
of the zone boundary magnon softening, of SSM045 
with that of Pr .63Sr .37MnO 3 (PSM037) @ and 
Lai-zSr^MnOa (x = 0.2,0.3; LSMO20, LSMO30) 0, 
the authors concluded that J4/ J\ increases dramatically 
for Ai_a;74.^Mn03 with x > 0.3. Since theoretical analy- 
sis based on the local density approximation + Hubbard 
U band calculations reveal that this doping dependence 
is consistent with the effect of rod-like (3z 2 — r 2 ) orbital 
correlations, the authors argue that the observed zone 
boundary magnon softening in Ai^ x A' x yhiOy, is due to 
the (3z 2 — r 2 )-type orbital fluctuations, in sharp contrast 
to all previous proposals [l5| . 

In this Letter, we take an approach different from that 
of Endoh et al., as we know that J4/J1 is nonzero for 
x = 0.3 manganites such as LSMO30 (Figs. 3 and 4) 
and La .7Cao.3Mn0 3 (LCMO30) @, in contrast with the 
expectation of Refs. 0,^1 ■ We decided to systematically 
analyze all existing spin-wave excitation data and take 
additional data in the FM metallic state of Ai- x A x MnOs 
at judicially selected x. We find that the low-q spin-wave 
stiffness D is insensitive to x while spin- wave excitations 
are systematically renormalized near the zone boundary 
with J4/J1 proportional to x. We also find on-site dis- 
order and lattice distortion that control Tc's of different 
x = 0.3 manganites have no effect on J4/J1, different 
from the expectations of the disorder effect theory . 
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FIG. 1: (a) Crystal structure of Ai- x A' x Mn03 with mag- 
netic exchange coupling indicated, (b) Phase diagram of 
Pro.7Cao.3Mn03 in the T-H plane 19]. (c) Phase diagram 
of Pro.55(Cao.85Sro.i5)o.45Mn03 in the T-H plane from trans- 
port measurements. Our neutron experimental conditions are 
marked as the (red) upper triangle and the (blue) square in 
the phase diagrams. 



These observations cannot be explained consistently by 
any current theory, thus suggesting that more than one 
mechanism is at play in determining the spin dynamical 
properties of the Ai-xA^MnOz manganites. 

For this study, we used single crystals of 
La .75Cao.25Mn0 3 (LCM025), Pr .7Cao.3Mn0 3 
(PCMO30) and Pro.55(Cao.85Sro.i 5 )o.45Mn0 3 (PC- 
SM045) grown by the traveling solvent floating zone 
technique. We chose these samples because they 
represent a large span in carrier concentrations. 
While LCM025 has a FM metallic ground state with 
Tn = 191 K 0, PCMO30 [H E3 and PCSM045 
|2 1| exhibit AF insulating behavior at zero field but 
can be transformed into FM metallic phases by field 
cooling from room temperature (Fig. 1). Our neutron 
scattering experiments were performed on triple-axis 
spectrometers at the High-Flux Isotope Reactor (HFIR), 
Oak Ridge National Laboratory and the NIST Center 
for Neutron Research (NCNR). The momentum trans- 
fers q = (qx,Qy,Qz) in units of A -1 are at positions 
(h,k,l) — (q x a/2ir, q y b/2Tr, q z c/2ir) in reciprocal lattice 
units (rlu), where lattice parameters for LCM025, 
PCMO30, and PCSM045 are a « b w c w 3.87 A, 
3.856 A, and 3.834 A, respectively. 

Figure 2 shows constant-^ scans at representative 
wavevectors for LCM025, PCMO30 and PCSM045 
along the [1,0,0] direction. The excitation peaks are 
sharp and resolution-limited at low-q, but become weak 
in intensity near the zone boundary. To obtain the 
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FIG. 2: The (/-dependent spin-wave excitations in LCM025 
(a), PCMO30 (b), and PCSM045 (c) at low temperatures. 
The data at different q's are incrementally shifted for clarity. 
The horizontal bars are the instrumental resolution and small 
shoulders around 25 meV in (c) are phonon scattering. E 
versus q 2 are plotted in the insets to determine D. 



strength of the average magnetic interaction, we ana- 
lyze the low-g data using E = A + Dq 2 , where E is the 
spin- wave energy obtained from Gaussian fits, A is field- 
induced Zeeman gap [2(| and D is the spin-wave stiffness. 
The slopes of the E versus q 2 lines, shown in the insets of 
Fig. 2, yield D values of 150 ± 3, 145 ± 8, 152 ± 3 meVA 2 
for LCM025, PCMO30 and PCSM045, respectively. It 
is remarkable that all three samples exhibit very similar 
low-g behavior after they are driven into FM states either 
by temperature or by magnetic field. This suggests that 
the average kinetic energy derived from the hopping of 
the itinerant electrons between adjacent manganese ions 
is independent of carrier concentration. 

To determine the evolution of magnetic excitations in 
Ai^ x A' x yiii03 as a function of x, one must first under- 
stand the effect of on-site disorder arising from the mis- 
match between rare- and alkaline-earth ions, as such dis- 
order might induce anomalous spin dynamical behavior 



[13j . The disorder is characterized by a 2 = 



-f 2 ), 



where Xi is the fractional occupancies of A-site species, 
ri and r = JT xm are individual and averaged ionic 
radius respectively |22j,|23|. Fig. 3(a) summarizes spin- 
wave dispersions along the [1,0,0] direction for a se- 
ries of Ai_ x ^Mn0 3 with x « 0.30 0, i, ©, H 0, 
while the doping dependence of magnon excitations is 
shown in Fig. 3(b) [a ^3- The solid lines in the figure 
are phenomenological fits to the data using Heisenberg 
Hamiltonian E(q) = A + 25[J(0) - J2j Jije^*-^} 
with nearest neighbor (Ji) and fourth-nearest neighbor 
( J4) exchange coupling. In the low-g limit, E(q) = A + 
87r 2 5(Ji + AJ^)q 2 . This simple Hamiltonian gives a sat- 
isfactory description of the data, where J4/J1 measures 
the magnitude of the effective zone boundary magnon 
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FIG. 3: (a) Spin-wave dispersion curves of various 
Ao.yA'q 3 Mn03 manganites along the [£, 0, 0] direction, (b) 
Dispersion curves for a series of Ai- x A' x MnOi as a function 
of x. The solid lines are least-square fits using the Heisenberg 
model with Ji and J4. 



softening [lq . We note that our previous experience |g 
showed that J2 and J3 have small contributions to the 
total magnon dispersion. 

We are now in a position to determine the effect of 
on-site disorder (cr 2 ) and average ionic size (f) on spin- 
wave excitations of Aq^A'q 3 MnC>3. Figures 4(a-c) show 
the cr 2 dependence of the spin- wave stiffness D, the near- 
est neighbor exchange coupling J\ and the ratio of J4/ ' J\ 
in Ao.jAq 3 MnC>3. With increasing disorder, the long- 
wavelength limit spin- wave stiffness shows no systematic 
trend but falls within a bandwidth of D = 160 ± 15 
meVA 2 [Fig. 4(a)]. While such behavior at low-c/ is not 
unexpected |7|, numerical calculations suggest a signif- 
icant zone boundary magnon softening with increasing 
disorder [l3| . In other words, increasing cr 2 should have 
no effect on D but dramatically increase J4/J1. Surpris- 
ingly, Figure 3(a) reveals no direct correlation between 
cr 2 and zone boundary magnon energy; and Figure 4(c) 
shows that J4/J1 is independent of cr 2 . Therefore, the 
on-site disorder has no observable effect on zone bound- 
ary magnon softening. 

In addition to inducing on-site disorder, replacing A 
by A' in A$/jA' 3 MnC>3 will also change the average ionic 
radius f and modify the length and angle of Mn-O-Mn 
bonds, thus leading to changes in effective transfer inte- 
grals between Mn ions or the bandwidth of the electrons 
|l7| . Figures 4(d-f) show the f dependence of the spin- 
wave stiffness D, 2SJi, and J4/J1. With increasing f, 
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FIG. 4: The on-site disorder dependence of (a) the spin-wave 
stiffness D obtained from low-g excitations, (b) the coupling 
25 Ji and (c) the ratio of Ja/Ji- The average ionic radius 
dependence of (d) D, (e) 2SJi, and (f) J4/J1. Dashed lines 
are guides to the eye. 



D shows a parabolic curve within a small bandwidth. It 
increases from 160 meVA 2 at f - 1.21 A for NSMO30 to 
176 meVA 2 at 1.25 A for LSMO30, and then decreases 
to 152 meVA 2 at 1.29 A for LBMO30 [Fig. 4(d)]. The 
Tc 's for NSMO30, LSMOSO^and LBMO30 are 198 K 
0, 351 K 0, and 350 K 40] respectively. Although 
increasing f leads to rapid changes in Tc, the kinetic 
energy (D) or bandwidth of the electrons only changes 
slightly Furthermore, J4/J1 is independent of f 

[Fig. 4(f)], thus indicating that the magnitude of the zone 
boundary magnon softening is independent of Tc and a 
general feature of the Aq^A' q 3 Mn03 manganites. 

Assuming the effect of on-site disorder and ionic size is 
weakly doping dependent, we can then study how spin- 
wave excitations of Ai- x A' x Mn03 are modified as a func- 
tion of x. For LSMO20, we used the stiffness value of 
D = 166.8 ± 1.51 meVA 2 obtained by high-resolution 
cold neutron triple-axis at low-g ^(| because this value 
is more accurate than the earlier result of 120 meVA 2 ob- 
tained on a thermal triple-axis || . Figure 5 summarizes 
the x dependence of D, 2SJ\, and J4/J1. Surprisingly, 
the spin- wave stiffness D is around 160 ± 15 meVA 2 and 
essentially unchanged for 0.2 < x < 0.45 [Fig. 5(a)], 
while T c varies from 305 K for LSMO20 16] to 135 K 
for SSM045 [l5|- This is in sharp contrast to the expec- 
tation of all DE based models, where D increases with x 
[solid line in Fig. 5(a)] 0- This also differs from a con- 
ventional ferromagnet, where Tc should be proportional 
to D. On the other hand, 2SJ\ decreases and J4/ Ji in- 
creases, approximately linearly, with increasing x [dashed 
lines in Figs. 5(b,c)]. Therefore, J4/J1 does not exhibit 
a huge rise in magnitude for x > 0.4 as expected from 
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FIG. 5: Doping dependence of (a) spin-wave stiffness D, (b) 
nearest neighbor exchange coupling Ji and (c) ratio of J4/J1. 
Dashed lines are guide to the eye. The solid line in (a) is the 
prediction of Ref. [3]. The (green) solid, (blue) dash-dotted, 
and (red) dashed lines in (c) are calculations from different 
models in Ha. 



the (3z 2 — r 2 )-type orbital fluctuations (dashed line), nor 
does it follow the predictions of the (x 2 — y 2 )-type orbital 
fluctuations (dash-dotted line) or free hybridized bands 
(solid line) shown in Fig. 5(c) flj) . 

Our systematic investigations in Figs. 4 and 5 put 
stringent constraints on microscopic theories of the zone 
boundary magnon softening. The possibility of on-site 
disorder- induced zone boundary softening is ruled out, 
as such a theory expects an enhanced softening with ei- 
ther increasing disorder or decreasing x |l3| . both con- 
trary to the observation. Similarly, it is unclear how on- 
site Coulomb repulsion in a DE mechanism can explain 
the doping independent behavior of the spin-wave stiff- 
ness Q- The magnon and Ai- x A' x -s\te optical phonon 
coupling (crossing) scenario postulated for LCMO30 8] 
also has difficulty in explaining the evolution of J4/J1, as 
the average A\- x A' x -aite mass and frequencies of associ- 
ated optical phonon modes do not vary dramatically from 
LCM025 to SSM045. If the large J4/J1 for x = 0.45 
materials stems from fluctuations of the (3z 2 — r 2 ) or- 
bital, J4/J1 should have a spectacular doping dependence 
around x = 0.4 [TH)]. However, this is not observed. Fur- 
thermore, the (3z 2 — r 2 )-type orbital fluctuations should 
have little or no effect on spin-wave softening at low 
carrier doping of x = 0.25 and 0.3. Therefore, it can- 
not be the origin of zone boundary magnon softening in 
Ai_ x A' x Mn03 at all doping levels. Finally, we note that 
the free hybridized band model and (x 2 — y 2 ) orbital fluc- 
tuation effects do not have the correct doping dependence 
for J4/J1 [Fig. 5(c)]. Although none of the current theo- 



ries is capable of explaining the spin dynamics in the en- 
tire doping regime, it is possible that more than one effect 
determines the properties of spin excitations in mangan- 
ites. For example, strong Coulomb repulsion in connec- 
tion with the electron-phonon coupling can induce orbital 
polarization throughout the phase diagram, ft is also 
known that chemical disorder can fundamentally modify 
the ground states of manganites an d be responsible 
for changes in spin dynamics . To understand the evo- 
lution of spin dynamical behavior in Ai^ x A' x yh\0^, one 
must consider interactions among spin, charge, orbital, 
and lattice degrees of freedom. 
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